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Abstract 

In this work, FRW universe filled with dark matter (perfect fluid with negligible pressure) along 
with dark energy in the background of Galileon gravity is considered. Four dark energy models with 
different EoS parametrizations have been employed namely, Linear, CPL, JBP and Logarithmic 
parametrizations. From Stern, Stern+BAO and Stern+BAO+CMB joint data analysis, we have 
obtained the bounds of the arbitrary parameters wo and ui 1 by minimizing the y 1 2 3 test. The 
best-fit values and bounds of the parameters are obtained at 66%, 90% and 99% confidence levels 
which are shown by closed confidence contours in the figures. For the logarithmic model unbounded 
confidence contours are obtained and hence the model parameters could not be finitely constrained. 
The distance modulus fi( z) against redshift z has also been plotted for our predicted theoretical 
models for the best fit values of the parameters and compared with the observed Union2 data sample 
and SNe Type la 292 data and we have shown that our predicted theoretical models permits the 
observational data sets. From the data fitting it is seen that at lower redshifts (z < 0.3) the 
SNe Type la 292 data gives a better fit with our theoretical models compared to the Union2 data 
sample. So, from the data analysis, SNe Type la 292 data is the more favoured data sample over 
its counterpart given the present choice of free parameters. From the study, it is also seen that 
the logarithmic parametrization model is less supported by the observational data. Finally we 
have generated the plot for the deceleration parameter against the redshift parameter for all the 
theoretical models and compared the results with the work of Farooq et al, 2013. 
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1 Introduction 


Over the past decade the most pressing goal for modern cosmology has been to find a proper explanation for 
the recent cosmic acceleration as confirmed from the observations of type la Supernovae and Cosmic Microwave 
Background (CMB) [H ©US SI 0 radiation. The standard explanation invokes an unknown component which has 
the property of positive energy density and negative pressure, known as “dark energy”. Different observations 
have also been shown that the contributions of dark energy and dark matter are respectively about 70% and 26% 
of the total energy of the universe. The phenomenon of the accelerated expansion of the universe has also been 
strongly confirmed by some other independent experiments like Sloan Digital Sky Survey (SDSS) [6], Baryonic 
Acoustic Oscillation (BAO) 0, WMAP data analysis [9, 10J etc. Over the past decade, there have been many 
theoretical models which mimic the dark energy behaviors. The simplest of them is the cosmological constant 
in which the equation of state is independent of the cosmic time and which fit the observations well. This model 
is the so-called ACDM model, containing a mixture of cosmological constant A and cold dark matter (CDM). 
However, this model is plagued by two famous problems, namely “fine-tuning” and the “cosmic coincidence” 
problem. In order to solve these two problems, numerous dynamical dark energy models were suggested, whose 
equation of state evolves with cosmic time. The scalar field or quintessence mm is one of the most favourite 
candidate of dark energy which produce sufficient negative pressure to drive the cosmic acceleration. In order 
to alleviate the cosmological-constant problems and explain the accelerated expansion, many dynamical dark 
energy models have been proposed, such as K-essence, Tachyon, Phantom, quintom, Chaplygin gas model, etc 
da nn mum in]. Interaction between dark energy and dark matter have been identified as a very important 
tool to address these problems. As a result interacting dark energy models including Modified Chaplygin gas 
[T5] , holographic dark energy model m and braneworld model |20] have been proposed. Since the discovery of 
General theory of relativity, the most pressing objective for cosmology has been found as the complete theory 
of quantum gravity i.e., the theory of everything. Relentless efforts and their results can be found widely in 
literature. Recently, such efforts using the principle of quantum gravity, led to the proposals of agegraphic dark 
energy (ADE) and the new agegraphic dark energy (NADE) models by Cai [21] and Wei et al. [22] respectively. 
Many theoretical models have been tallied with the observations using different data sets, say TONRY, Gold 
sample data sets [221LS (HI US] etc- In Einstein’s gravity, the modified Chaplygin gas [TS] best fits with the 
3 year WMAP and the SDSS data with the choice of parameters A = 0.085 and a = 1.724 [26] which are 
improvements over the previous ones —0.35 < A < 0.025 [27]. 

Another possibility of the accelerated expansion is attributed to the fact that general relativity is only 
accurate on small scales and therefore needs modification on cosmological distances. This led to the proposal of 
modified gravity theories. In this case, cosmic acceleration would arise not to from dark energy as a constituent 
substance but rather from the dynamics of modified gravity. Modified gravity constitutes an interesting dy¬ 
namical alternative to ACDM cosmology regarding given the fact that it is also able to describe the current 
acceleration in the expansion of our universe. One of the simplest modified gravity is DGP brane-world model 
[251 . The other alternative approach dealing with the acceleration problem of the Universe is changing the 
gravity law through the modification of action of gravity by means of using f(R) gravity [29jj30] instead of the 
Einstein-Hilbert action. Some of these models, such as 1 /R and logarithmic models, produced an accelerated 
expansion for the Universe at the present time [3Tj. Other modified gravity includes /(T) gravity, /(G) gravity, 
Gauss-Bonnet gravity, Horava-Lifshitz gravity, Brans-Dicke gravity, etc [32] [331 SU GUI 36. • 

Of late an infrared modification of classical gravitation was proposed, which is a generalization of the 4D 
effective theory in the DGP model m- The theory considers a self-interaction term of the form Cl <f> (Vt/>) 2 in 
order to recover GR in high density regimes. The most notable feature of the theory is that it is invariant under 
the Galileon shift symmetry S^(f> —> S^cf) + in the Minkowski background. Because of this invariance, the 
equation of motion remains a second order differential equation, preventing the introduction of extra degrees of 
freedom, which are usually associated with instabilities. So we assume the FRW universe in Galileon gravity 
model filled with the dark matter and dark energy. 

As stated earlier, in order to explain the evolution of the universe, various dark energy models have been pro¬ 
posed. All these models must be constrained by astronomical observations in order to bring about a consistency 
between the proposed theoretical models and the observational data. For any dark energy model, the Equation 
of state (EoS) parameter ui plays a vital role and can reveal the nature of the model which is responsible for the 
cosmic acceleration. Different form of EoS lead to different dynamical nature and may influence the evolution 
of the Universe considerably. The EoS parameter to and its time derivative with respect to Hubble time are 
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currently constrained by the distance measurements of the type la supernova. From literature it is seen that the 
current observational data constrain the range of EoS as —1.38 < ui < —0.82 [38]. Recently, the combination 
of WMAP3 and Supernova Legacy Survey data showed a significant constraint on the EOS to = — 0.97^09 for 
the dark energy in a flat universe 139] . Recently, some parametrizations for the variation of EOS parameters 
ui(z) have been proposed describing the DE component. We will discuss them in detail in section 3. 

It is a known fact that the standard big bang model of cosmology is crippled by various problems such as 
the horizon problem, the flatness problem, the singularity problem, etc. Considering the inflationary model, 
we have been able to explain the flatness problem. Now for a flat universe filled with matter and dark energy, 
an accurate knowledge of f2<j e and u>d e for dark energy can only be obtained if we have a proper knowledge of 
Q rn for dark matter and H(z) [231 SO] . For z > 0.01, we see that the TONRY data with 230 data points [24] 
with 23 data points from Barris et al. [25] is quite acceptable. In the redshift range 1 < z < 1.6, the GOLD 
sample of Riess et al. [3] with 156 data points is valid. From the CMBR data one finds that Da + = 1 TO] , 

From the most recently obtained Riess data set the best fit value of fl m is obtained as 0.31 ± 0.04. This value 
matches with the value of Q m = 0.291q'q 3 obtained from the previous Riess data [2]. In Chowdhury et al. [25] 
the value was found to be 0.31 ± 0.08. For a flat universe, the best fit value for the equation of state parameter 
w for union 2 data sample is w = —0.997 ±o:o54( stat )to ;gg 2 (stat+sys together). For a curved one, the value is 
w = —1.038ig o5g(stat)to;o9r( sta t+ s y s together) [4lj . 

The success of any dark energy or modified gravity model, depends basically on its consistency with the 
observational data. This is our basic motivation for the work. We reconstruct the hubble parameter H using 
the parameters of dark energy, dark matter and modified gravity. Then we set up a comparison scenario 
between the reconstructed H (Htheoretical) and the values of H obtained from observational data (H 0 b serva u 0 na i). 
This is accomplished by the procedure of chi-square test. Recently, numerous works can have appeared in 
the literature, which aims at constraining dark energy model parameters using Stern and union2 data sets 
@21031 Banana- The basic concepts of Galileon gravity theory are presented in section 2. The behaviour of 
some reconstructed cosmological parameters is studied in section 3. Some basic calculations using the Galileon 
gravity and parametrized dark energy models is given in section 4. The observational data analysis tools in 
observed Hubble data (OHD) or H(z)-z (Stern), OHD+BAO and OHD+BAO+CMB for y 2 minimum test will 
be studied in section 5 and wc will also investigate the bounds of unknown parameters (wo,wi) for the various 
dark energy parametrizations by fixing other parameters. The best-fit values of the parameters are obtained at 
66%, 90% and 99% confidence levels. The redshift magnitude observations for our theoretical parametrization 
models in Galileon gravity for the best fit values of the parameters and the observed SNe la union2 data sample 
is studied in section 6. Analysis of our theoretical models with SNe Type la 292 data is given in section 7. 
Finally the paper ends with a discussion in section 8. 


2 Basic equations of Galileon gravity 


The Galileon gravity theory is described by the action [371 HZ1 mi HU [50]: 


S = 


d A Xy/~^g 


<j)R 


~r (V</>) 2 + /(</>)□ 0 (V</>) 2 + £ m 
<P 


(1) 


where (f> is the Galileon field, w is known as Brans-Dicke (BD) parameter, the coupling function /(</>) has 
dimension of length, (Vqf>) 2 = v (f), □ </> = and C m is the matter Lagrangian. Variation of 

Eqn. (1) with respect to the metric g gives the Einstein’s equations, 


w 


G- = # + 4 , (v ^ “ ^ U(t))+ ^ 

l ( l 


V„(/) - (V0) 2 


- - 2(V^) 2 ]VV - V„[/(</>) + /(^)V^V^ □ <j> 


( 2 ) 


For the Friedmann-Robertson-Walker background metric, the Einstein’s field eqns for Galileon gravity gives, 


^n 2 = {-- + f - 12 + t 2 fW (zh - ^/) i 3 


(3) 
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- 3H 2 - 2H = + J + I 2 + 2 HI + -I 2 - </ 2 /(</>) (i + 2 + ai I 2 ] 1 2 (4) 

2(p 2 \ 2 / 

where H(t) = %, I(t) = $ and a„ [<£(*)] = . 


Here p = + p m and p = +Pm, where p m , p m are the energy density and pressure of the dark matter 

and pde, Pde are respectively the energy density and pressure contribution of some dark energy. 


3 Some parametrizations on EoS parameter 

Recently some parametrizations for the variation of EOS parameters uj(z) have been proposed describing the 
DE component. Here we discuss some of them. 

3.1 Linear Parametrization 

Here the EoS is given by EB 


ui(z) = wo + u>iz (5) 

Here wo = —1/3 and wi = —0.9 with z < 1 when Einstein gravity has been considered. This grows increasingly 
unsuitable for 2 > 1. 


3.2 Chevallier-Polarski-Lindler (CPL) Parametrization 


Here the EoS is given by 


uj(z) = w 0 + Wi 


z 

1 + z 


( 6 ) 


This ansatz was first discussed by Clrevallier and Polarski [52| and later studied more elaborately by Linder 
[53] . In Einstein gravity the best fit values for this model while fitting with SNla gold data set are wo = —1.58 
and wi = 3.29. This parametrization has several advantages over its predecessor. Some of them are: 1) It has 
a manageable 2-dimensional phase space; 2) reduces to the old linear redshift behaviour at low redshift; 3) has 
a well behaved and bounded behaviour for high redshift; 4) highly accurate in reconstructing many scalar field 
equations of state and the resulting distance-redshift relations; 5) highly sensitive to observational data and 
6) simple physical interpretation. In addition to the above advantages, the new parametrization is also more 
accurate than the previous one. The most remarkable thing about the parametrization is that, it reconstructs 
the distance-redshift behaviour of the SUGRA model [54] to 0.2% over the entire range out to the last scattering 
surface (z ~ 1100). 


3.3 Jassal-Bagla-Padmanabhan (JBP) Parametrization 

Here the EoS is given by [55] 

w(z) =W 0 +W l ( 1+z )2 ( 7 ) 

A fairly rapid evolution of this EoS is allowed so that w(z) > —1/2 at z > 0.5 is consistent with the supernovae 
observation in Einstein gravity. From the EoS it can be seen that w(0) = wo and w(oo) = wo- For a standard 
cosmological model with a hot big bang to be valid, we must have u(z 1) < —1/3. 


3.4 Logarithmic Parametrization 

Here the EoS is given by EH1E7] 

w(z) = w 0 + Wi log(l + z) (8) 

This model is sparsely found in literature, and therefore very little is known about the behaviour of this model. 
It has been used now and then for parametrization of dark energy. Here we will analyze it and try to compare 
it with other models. This will help us to understand its importance in cosmology and also to find out the 
reason why it has been so seldom studied in past. This is the basic reason for the inclusion of this model in the 
present study. 
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4 Basic calculations on Galileon gravity using various EoS parametriza- 
tions 

Here we will consider the flat, homogeneous and isotropic universe described by the FRW background metric. 

The corresponding Einstein equations in Galileon gravity are given by eqns.(3) and (4). The DE equation of 
state is taken as pd e = <x(z)pd e - We neglect the pressure contribution of matter, i.e. p m = 0. The independent 
conservation equations for DE and DM are given as, 



pde + 3 H (Pde + Pde) — 0 

(9) 

and 


Pm + 3 H Pm — 0 

(10) 

From eqn.(10), we have 


Pm — PmO (1 Z) 

(11) 


where we have used the relation a = z being the cosmological redshift parameter. Here p m o is a constant. 
Integrating equation (9) for the four DE models we get four different solutions for DE density. Now, we consider 
dimensionless density parameters D m0 = and fldeO = and for simplicity choosing f{cj)) = fo<j> n and 
(j) = <j)oa m (fo > 0, (f>a > 0, m > 0, n > 0) in the power law form, [43, we have the expression for Hubble 
parameter H in terms of redshift parameter z as follows (from eq. 3): 


H(z) 


1 

2 (6f 0 m 3 4>l +n - fom 4 rut>Q +n ) 



6rn-6)(l + z) m(2+n Vo 


([m 2 w - 6m - 6)(1 + z) m ^ 2+n Uo) 2 - / 0 mVg +n - 12H 2 ((1 + z ) 3+m ^n mo + 




( 12 ) 


4.1 Expressions for Linear parameterization 

Using equation (5), we integrate equation (9) and get the expression for DE density as 

Pde = Pde o (1 + exp |3 Wl (i + z'jj (13) 

where pd e o is a constant. In terms of the above mention dimensionless density parameters and redshift parameter 
(z) we calculate the Hubble parameter as follows, 


H-lin (z) — 


3 (1 + z ) m ( 2 +-)+ 3 -i fo (i + m _ y w ) 
fom 3 (/)Q +n (1 + z) 3ui (6 — mn) 


^36 (1 + .) 2 ™(2+")+6-1 $ (l + m 2 - im 2 ^" 


+12H; 


((1 + 0 ) m+m(2+ ") +3 “ 1 n m0 (1 + 3^ + 3z 2 + z 3 ) + exp{3u! (1 + 0 )} (1 + z )™+-( 2 +G+3- 0 ^ 


(l + 3z + 3 0 2 + 0 3 )) f Q m 3 (1 + z) 3uJl 


<t>l +n 


(■mn — 6)^j 



2/ 0 m 3 (1 + z) 3ui1 (j) Q +ra (mn — 6) 


(14) 
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4.2 Expressions for CPL parameterization 

For this DE parametrization the DE density is obtained as, 


Pde = Pde 0 (1 + zf {1+U>0+Ull) exp { 3Wl (1 + 2:)} 


(15) 


Using the dimensionless density parameters, the expression for the Hubble parameter in terms of the redshift 
is obtained as, 


Hcpl(z) = 


3 (1 + 2) I7 b 2+J 4 <^ 0 (l + jTi — 1 m 2 w) 
/om 3 (/>g +rl (6 — mn) 


2/ o m 3 0 o + " (6 — mn) 


36 (1 + z) 


2m(2+n) 


1 + m - m z w 

6 


+ 121 + 


; ((1 + 2:) m+m(2+n) fimO (1 + 32 + 32 2 + 2 3 ) + CXp | | ( X + Z ) 

1/2 X 

(l + 32 + 32 2 + 2 3 )) f 0 m 3 (1 + z) 3ui1 (mn - 6)^ j 


3^1 \m+m(2+n)+3a;o+3a;i q 


de 0 


1 1/2 


2/o?n 3 </>o + " (mn — 6) 


(16) 


4.3 Expressions for JBP parametrization 

Here the DE density is obtained as, 


Pde = Pde o (1 + 2) 1 Uo) exp < - 


3 (l+«o)_/ 3aJl ( g +l) 


(1 + 2) 2 


(17) 


Using the dimensionless density parameters, the expression for the Hubble parameter in terms of the redshift 
is obtained as, 


Hjs p (z) [\^/ oTO 3(_6 + mn) 


exp 


3(^1 + 2zUi) \ ,-3-r 

2(1 + 2) 2 J 00 


exp 


3(o/i -|- 2,z(jj 
2(1 + z)“ 


— | (—6 + m (—6 + mw )) 


(1 + 2) m(2+n Vo ( exp 


3(o/i + 2iZuj±^ 

(1 + ^) 2 


3(c»;^ -|-22 :cj j ) 

(6 + m(6 — mw)) 2 (l + 2) 2m ^ 2+ "Vo+ 12e 2 < 1+2 + /oHo?ti 3 (—6 + mn) 


(1 + *) 3 +"d 3 +^ 3 +« (exp j P m0 + (1 + z) 3 “°n x0 


1/2N 


1 / 2 ' 


(18) 


4.4 Expressions for Logarithmic parametrization 

For this parametrization model, the DE density is obtained as, 

Pde = PdeO (1 + *) 3 ( 1+ -°+ Wl ) (19) 

Using the dimensionless density parameters, the expression for the Hubble parameter in terms of the redshift 
is obtained as, 


HLog{z) — 


LV2/om 3 (—6 + mn) 


4 3 " (—6(1 + 2)”^ 2+ "Vo — 6m(l + 2) m ^ 2+ "Vo + m 2 w (1 + z ) m( ' 2+n ' 1 + ((6 + 6 m 


—m 2 w) 2 (1 + 2) 2m(2+n Vo + 12/offo 2 ™ 3 (-6 + mn) (1 + z ) 3+m ^ ^ (p m0 + (i + ^) 3 (-o+-i) 


1/2 


1 1/2 


( 20 ) 
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5 Observational data analysis mechanism 

In the following subsections, we present the observational data analysis mechanism for Stern, Stern+BAO and 
Stern+BAO+CMB observations. We use the y 2 minimum test from theoretical Hubble parameter with the 
observed data set and find the best fit values of unknown parameters for different confidence levels. From the 
above expressions, we see that H(z ) contains the unknown parameters like u>o, Wi, f2 m o, fi deO , w, n, m, /o, 
</>o- Now the relation between two parameters will be obtained by fixing the other parameters and by using 
observational data set. Eventually the bounds of the parameters will be obtained by using this observational 
data analysis mechanism. 


z 

H{z) 

a(z) 

0 

73 

± 8 

0.1 

69 

± 12 

0.17 

83 

± 8 

0.27 

77 

± 14 

0.4 

95 

± 17.4 

0.48 

90 

± 60 

0.88 

97 

± 40.4 

0.9 

117 

± 23 

1.3 

168 

± 17.4 

1.43 

177 

± 18.2 

1.53 

140 

± 14 

1.75 

202 

± 40.4 


Table 1: The Hubble parameter H(z) and the standard error a(z) for different values of redshift 2 . 

5.1 Analysis with Stern data 

We analyze the model, using observed value of Hubble parameter at different redshifts (twelve data points) 
listed in observed Hubble data by Stern et al. [58]. The Hubble parameter H(z) and the standard error er(z) for 
different values of redshift z are given in Table 1. Since we are using testing of hypothesis, so before proceeding, 
we form our null and alternate hypothesis which are given below. 


Null Hypothesis: Hqi Hth eore tica.i — Hobservational 

Alternative Hypothesis: Hibeoreticai ^ Hobservational 


Here we test the null hypothesis Hq against the alternative hypothesis H\. For this purpose we first form 
the \ 2 statistics as a sum of standard normal distribution as follows: 


\Stern 


E 


(H(z) - H obs (z)) 2 
a 2 (z) 


( 21 ) 


where H(z) and H a b s (z ) are theoretical and observational values of Hubble parameter at different redshifts 
respectively. Here, H 0 b s is a nuisance parameter and can be safely marginalized. We consider the present value 
of Hubble parameter Hq = 72 ± 8 Kms -1 Mpc -1 and a fixed prior distribution. Here we shall determine best 
fit value of the parameters (u>i,uiq) by minimizing the above distribution Xstern an d fixing the other unknown 
parameters with the help of Stern data. 

We now plot the graph for different confidence levels. Now our best fit analysis with Stern observational data 
support the theoretical range of the parameters. The 66% (solid, blue), 90% (dashed, red) and 99% (dashed, 
black) contours are plotted in figures la, 16, lc and Id for VL m0 = 0.28, SldeO = 0.72, n = 3.155, m = 2, u> = —25, 
/o = 0.9, <j >0 = 0.1. The best fit values of (wi,wo) are tabulated in Table 2. Here the values of n m o and Qdeo 
are motivated from the recent cosmological observational data. The value of ui is chosen in accordance with 
reference |7j- The rest of the parameters namely, /o, <f>o, m and n arises due to our particular choice function 
of the scalar field. Due to their secondary importance we consider them as free parameters and their values are 
fine tuned in such a way so that effective closed contours are obtained. 
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5.2 Joint analysis with Stern+BAO datasets 

The method of joint analysis, the Baryon Acoustic Oscillation (BAO) peak parameter value has been proposed 
by 0 and we shall use their approach. Sloan Digital Sky Survey (SDSS) survey is one of the first redshift 
survey by which the BAO signal has been directly detected at a scale ~ 100 MPc. The said analysis is 
actually the combination of angular diameter distance and Hubble parameter at that redshift. This analysis is 
independent of the measurement of H$ and not containing any particular dark energy. Here we examine the 
parameters ui\ and u>o for parametrized DE model from the measurements of the BAO peak for low redshift 
(with range 0 < z < 0.35) using standard y 2 analysis. The error is corresponding to the standard deviation, 
where we consider Gaussian distribution. Low-redsliift distance measurements is lightly dependent on different 
cosmological parameters, the equation of state of dark energy and have the ability to measure the Hubble 
constant H$ directly. The BAO peak parameter may be defined by 


A = 


A I? 


E(z i) 1 / 3 


1 

Zl 


dz 

W) 


2/3 


( 22 ) 


Here E(z) = H(z)/Hq is the normalized Hubble parameter, the redshift zi = 0.35 is the typical redshift of the 
SDSS sample and the integration term is the dimensionless comoving distance to the to the redshift z\ The 
value of the parameter A for the flat model of the universe is given by A = 0.469 ± 0.017 using SDSS data 0 
from luminous red galaxies survey. Now the y 2 function for the BAO measurement can be written as 


Xbao ~ 


(A- 0.469) 2 
(0.017) 2 


Now the total joint data analysis (Stern+BAO) for the y 2 function may be defined by 


(23) 


Xtotal Xstern XbAO ( 2 ^) 

According to our analysis the joint scheme(Stern+BAO) gives the best fit values of {oji,ujq) in Table 3. Finally 
we draw the contours for the 66% (solid,blue), 90% (dashed, red) and 99%(dashed, black) confidence limits 
depicted in figures 2a, 2b, 2c and 2d for H m0 = 0.28, fldeO = 0.72, a = 0.001, n = 0.5,m = 10, oj = —3, 
w m = 0.03, /„ = 0.01, 0 O = 0.01. 


5.3 Joint analysis with Stern+BAO+CMB datasets 

One interesting geometrical probe of dark energy can be determined by the angular scale of the first acoustic 
peak through angular scale of the sound horizon at the surface of last scattering which is encoded in the 
CMB power spectrum Cosmic Microwave Background (CMB) shift parameter is defined by 05J 62) EH- It is 
not sensitive with respect to perturbations but are suitable to constrain model parameter. The CMB power 
spectrum first peak is the shift parameter which is given by 


n = Jn. 


dz 

W) 


(25) 


where 2 2 is the value of redshift at the last scattering surface. From WMAP7 data of the work of Komatsu et 
al. [62| the value of the parameter has obtained as 72. = 1.726 ± 0.018 at the redshift z = 1091.3. Now the y 2 
function for the CMB measurement can be written as 


2 (72 — 1.726) 2 

XcMB ~ (0.018) 2 


(26) 


Now when we consider three cosmological tests together, the total joint data analysis (Stern+BAO+CMB) for 
the y 2 function may be defined by 


XtOTAL — Xstern + X%AO + XCMB ( 27 ) 

Now the best fit values of (tt+wo) for joint analysis of BAO and CMB with Stern observational data support 
the theoretical range of the parameters given in Table 4. The 66% (solid, blue), 90% (dashed, red) and 99% 
(dashed, black) contours are plotted in figures 3a, 3 b, 3c and 3 d for S2 m0 = 0.28, fide o = 0.72, a = 0.001, 
n = 0.5 ,m = 10, ui = —3, w m = 0.03, fo = 0.01, <po = 0.01. 











Fig.la Fig.lb Fig.lc Fig.Id 

Figs la, lb, lc and Id show the variation of ujq with oj± for different confidence levels for linear, 

CPL, JBP and logarithmic models respectively. The 66% (solid, blue, the innermost contour), 90% 
(dashed, red, next to the innermost contour), and 99% (dashed, black, the outermost contour) contours 
are plotted in these figures for the H(z) — z (Stern) analysis. 



Fig.2a Fig.2b Fig.2c Fig.2d 

Figs 2a, 2b, 2c and 2d show the variation of ljq with oj\ for different confidence levels for linear, 


CPL, JBP and logarithmic models respectively. The 66% (solid, blue, the innermost contour), 90% 
(dashed, red, next to the innermost contour), and 99% (dashed, black, the outermost contour) contours 
are plotted in these figures for the H{z)z (Stern+BAO) joint analysis. 


DE models 

Best fit values 




w 0 

Wl 

Xmin 

Linear 

-1.96664 

0.992802 

11.4933 

CPL 

-1.65925 

0.967584 

8.4682 

JBP 

-0.932444 

-1.82202 

7.88769 

Log 

-0.716018 

0.5 

127.572 


Table 2: STERN: The best fit values of ojq and wi and the minimum values of \ 2 for different redshift 
parametrization models. 


DE models 

Best fit values 




w 0 

Wl 

V 2 ■ 

A mm 

Linear 

-10 

2.42671 

928.664 

CPL 

-2.39571 

1.76133 

773.265 

JBP 

-0.980504 

-3.12768 

781.012 

Log 

2.75676 

-3 

903.023 
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Fig.3a Fig.3b Fig.3c Fig.3d 

Figs 3a, 3b, 3c and 3d show the variation of uq with for different confidence levels for linear, 


CPL, JBP and logarithmic models respectively. The 66% (solid, blue, the innermost contour), 90% 
(dashed, red, next to the innermost contour), and 99% (dashed, black, the outermost contour) contours 
are plotted in these figures for the H(z)z (Stern+BAO+CMB) joint analysis. 




ojO 




Fig.4a Fig.4b Fig.4c Fig.4d 

Figs 4 a > 4b, 4 C an d 4d show the variation of ujq with ui\ for different confidence levels for linear, 


CPL, JBP and logarithmic models respectively. The 66% (solid, blue, the innermost contour), 90% 
(dashed, red, next to the innermost contour), and 99% (dashed, black, the outermost contour) contours 
are plotted in these figures for the H(z)z of SNe Type la 292 data. 


Table 3: STERN+BAO: The best fit values of wo and wi and the minimum values of y 2 for different redshift 
parametrization models. 


DE models 

Best fit values 




w 0 

Wl 

Xmin 

Linear 

-3.01522 

1.53659 

9969.52 

CPL 

-3.70734 

0.750595 

11645.9 

JBP 

-2.9022 

-1.03615 

10400.9 

Log 

9.85277 

-10 

10105.5 


Table 4: STERN+BAO+CMB: The best fit values of wo and wi and the minimum values of y 2 for different 
redshift parametrization models. 


DE models 

Best fit values 




w 0 

Wi 

Xmin 

Linear 

-2.70613 

0.509612 

879522.25 

CPL 

-2.63487 

0.508819 

867545.51 

JBP 

-2.3388 

-1.00341 

847931.41 

Log 

-2.13811 

0.5 

1.32293 x 10 B 


Table 5: H(z)-z SNe Type la 292: The best fit values of wo and wi and the minimum values of y 2 for different 
redshift parametrization models. 
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Figs 5a and 5b show the variation of g(z) vs z for linear, CPL, JBP and logarithmic parametriza- 

tions (solid lines) and their comparison with observational data (blue dots). The dots denote the 
Union2 sample data in fig 5a and the SNe Type la 292 data in fig 5b. 


6 Redshift-magnitude observations from supernovae type la Union2 
sample (from Amanullah et al. 2010) 

The Supernova Type la experiments provided the main evidence for the existence of dark energy. Since 1995, 
two teams of High -z Supernova Search and the Supernova Cosmology Project have discovered several type la 
supernovas at the high redshifts DO [H El ESI- The observations directly measure the distance modulus of a 
Supernovae and its redshift z [64j :_65]. Now, take recent observational data, including SNe la which consists of 
557 data points and belongs to the Union2 sample [41]. From the observations, the luminosity distance d^lz ) 
determines the dark energy density and is defined by 

d L (z) = (l + z)H 0 J^ (28) 

and the distance modulus (distance between absolute and apparent luminosity of a distance object) for Super¬ 
novas is given by 


g(z) = 5 log 10 


d L {z)/H 0 
1 M Pc 


+ 25 


(29) 


The best fit of distance modulus as a function p(z) of redshift z for our theoretical models and the Supernova 
Type la Union2 sample are drawn in figure 5a for our best fit values of u>i and uiq for (Stern) + BAO + CMB 
joint analysis as H m0 = 0.28, Qdeo = 0.72, a = 0.001, n = 0.5,to = 10, w = —3, w m = 0.03, /o = 0.01, fo = 0.01. 
Similarly the best fit trajectories for our theoretical models are obtained and compared with SNe Type la 292 
data in fig.5b. From the curves, we see that the theoretical parametrized DE models in Galileon gravity is in 
agreement with the union2 sample data and the SNe Type la 292 data. 


7 Analysis with Supernovae Type la 292 data (from Riess et al.(2004, 
2007) and Astier et al.(2006)) 


Here we have analyzed our theoretical parametrized DE models in the light of the Supernovae Type la 292 data 
mmm- The compiled data can be found in ref m- Bounds for the model parameters (wo,Wi) have been 
obtained by constraining the theoretical DE models with the SN Type la 292 data, after fixing the cosmological 
parameters in their viable range. Given below is the y 2 statistic which has been used as a minimization technique 
for our data analysis, 


X(SNeTypeIa) 


= E 


(- H{z)~H obs {z)f 
a 2 {z) 


(30) 
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where H 0 b s (z), cr(z) are obtained from the SNe Type la 292 data table and H(z) is the calculated value of the 
Hubble parameter obtained from our theoretical model. The probability distribution function is given as, 

L = J e - 5X(SAr«T # p B Ja)p(tf 0 ) dH 0 (31) 

where the prior distribution function for Hq is P(Hq). The best fit values of the model parameters (wo,wl) 
obtained for different theoretical models by the y 2 minimization technique is presented in table 5 along with 
their respective minimum y 2 values. Contours, which are showing the numerical range of the parameters, are 
obtained for different confidence intervals and are given in the figs. 4a, 4b, 4c and 4d. 



-1.0 -0.5 0.0 0.5 1.0 


z 

Fig. 6 

Fig 6 shows the variation of deceleration parameter q vs redshift z for linear (Blue), CPL (Red), 
JBP (Black) and logarithmic (Green) parametrizations. 


8 Discussions and conclusions 

In this assignment, we have proposed the FRW universe filled with DM (perfect fluid with negligible pres¬ 
sure) along with dark energy in the background of Galileon gravity. For describing DE components, four 
parametrizations have been proposed for the variations of EOS parameter ui(z). They are Linear, CPL, JBP 
and Logarithmic parameterizations. We have presented the Hubble parameter H in terms of the observable 
parameters O m o, fide o, Ho with the redshift z and the other parameters like u>o and uq. We have chosen the 
observed values of H m0 = 0.28, H de0 = 0.72 and H 0 = 72 Kms~ 1 Mpc~ 1 . From STERN data set (12 points), 
we have obtained the best fit values of the arbitrary parameters w 0 and uq (Table 2) by minimizing the y 2 
test and by fixing the other parameters. Similarly, the best fit values of uiq and uq are presented in Table 3 
and Table 4 for STERN+BAO and STERN+BAO+CMB respectively. In Table 5, the best fit values of the 
parameters are given using the Supernovae Type la 292 dataset. Likelihood contours have been drawn for 
different confidence levels (66%, 90%, 99%) constraining the parameters (uo,ui). In the figs, la, lb, lc and Id 
contours have been drawn for STERN data. In figs. 2a, 2b, 2c and 2d contours are drawn for STERN+BAO 
joint data analysis. STERN+BAO+CMB joint data have been used to generate the plots 3a, 3b, 3c and 3d. 
Contours for Supernovae Type la 292 data can be found in figs. 4a, 4b, 4c and 4d. Finally in figs. 5a and 5b, 
using the best fit values for our different theoretical parametrization models, g vs redshift (z) curves are plotted 
and compared with Union compilation data and SNe Type la 292 data respectively. 

Figs, la to Id have been generated for Stern data. In fig. la, which is the plot for linear model, for 66% 
confidence interval u 0 and uq lies between —2.39 to —1.51 and 0.78 to 1.16 respectively. Similarly for 90% 
confidence limit, the ranges for uq are —2.52 to —1.58 and those for uq are 0.74 to 1.18. For 99% confidence 
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limits the corresponding ranges for u >o are —2.67 to —1.48 and for uq are 0.68 to 1.22. From fig. lb for CPL 
parametrization we see that for 66% confidence limits wo lies between —2.12 to —1.28 and aq lies between 0.76 
to 1.12. For 90% confidence limit ujq lies between —2.19 to —1.21 and uq lies between 0.72 and 1.14. Similarly 
for 99% limits the ranges are —2.34 to —1.10 for ujq and 0.66 to 1.18 for ui In fig. lc, for JBP model, the 
corresponding ranges for w 0 and cq are —1.35 to —0.65 and —2.20 to —1.45 respectively, in 66% confidence 
limits. For 90% confidence limits the ranges are —1.50 to —0.60 for ujq and —2.22 to —1.38 for cq. Finally for 
99% confidence interval the values are —1.65 to —0.55 for oj 0 and —2.23 to —1.27 for wi- 

In fig. 2a for linear model (joint analysis with Stern+BAO), the ranges of parameters are —11.50 to 0.50 
for w 0 and 0.12 to 2.70 for cq in 66% confidence interval. For 90% confidence level, w 0 varies from —12.30 to 
0.61 and ui i varies from —0.20 to 2.80. Finally for 99% confidence range, wo varies from —13.20 to 0.90 and cq 
varies from —0.25 to 2.85. In fig. 2b, for the CPL model, wo varies from —2.86 to —1.89 and cq varies from 1.47 
to 2.25 in 66% confidence interval. In 90%, confidence limits, uiq varies from —2.96 to —1.80 and cq varies from 
1.42 to 2.70. For 99% confidence limits the variation is from —3.08 to —1.68 for coo and 1.30 to 2.11 for cq. Plot 
for JBP is generated in fig. 2c. We see that in 66% confidence interval, loq varies from —1.31 to —0.69 and cq 
varies from —3.61 to —2.59. In 90% confidence limits, u>o varies from —1.39 to —0.61 and cq varies from —3.72 
to —2.48. Similarly in 99% confidence limits, ujq and cq varies respectively from —1.49 to —0.56 and —3.82 to 
-2.30. 

From figs. 3a to 3d, plots are generated for Stern+BAO+CMB joint data analysis. In fig. 3a, which is for 
linear model, the acceptable range for w 0 is obtained as —3.59 to —2.52 and that for cq is from 1.30 to 1.73 
at 66% confidence level. Similarly at 90% confidence level the values of u>o vary from —3.76 to —2.42 and for 
wi they vary from 1.26 to 1.76. Finally for 99% confidence interval, wo is from —3.96 to —2.32 and cq is from 
1.17 to 1.81. In fig. 3b, for CPL model, u o varies from —4.95 to —3.05 and cq varies from 0.61 to 0.81 in 66% 
confidence level. In 90% confidence limits the ranges of u >o and cq are respectively from —5.45 to —2.95 and 
0.55 to 0.85. Similarly for 99% limits, the acceptable range is from —6.05 to —2.82 for ojq and from 0.49 to 0.89 
for oj i. In fig. 3c, the plot has been generated for JBP model. For 66% confidence interval, ujq takes the values 
from —4.8 to —2.0 and wi takes the values from —1.4 to —0.45. At 90% confidence interval, the acceptable 
ranges for wo and cq are —5.4 to —1.9 and —1.50 to —0.43 respectively. Finally at 99% confidence limits wo 
varies from —6.5 to —1.8 and cq varies from —1.7 to —0.25. 

Plots for SNe Type la 292 data are generated in figs. 4a to 4d for different parametrization models. In fig. 
4a, for linear model, at 66% confidence level, loq varies from —2.7142 to —2.6980 and cq changes from 0.5078 
to 0.5116. At 90% confidence limits w>o ranges from —2.7159 to —2.6964 and cq ranges from 0.5074 to 0.5118. 
Finally at 99% confidence interval, the acceptable ranges are —2.7165 to —2.6945 for w 0 and 0.5070 to 0.5123 
for cq . Fig. 4b is generated for the CPL model. Here at 66% confidence interval, wo ranges from —2.6425 to 
—2.6270 and ui i ranges from 0.5070 to 0.5105. At 90% confidence limits, —2.6442 to —2.6247 is the range for 
wo and 0.5066 to 0.5111 is the range for cq. Similarly at 99% confidence interval, the acceptable ranges vary 
from —2.6460 to —2.6239 for wo and 0.5062 to 0.5116 for cq. Plot for the JBP model has been obtained in fig. 
4c. Here we see that, coq and cq respectively range from —2.3460 to —2.3318 and —1.0072 to —0.9999 for 66% 
confidence limit. At 90% confidence limits ojq attains values from —2.3475 to —2.3302. Similarly the acceptable 
range for cq is from —1.0078 to —0.9990. Likewise, at 99% confidence limits the ioq lies in the range —2.3490 
to —2.3287 and cq lies in the range —1.0088 to —0.9982. 

For all the four cases the Logarithmic parametrization is studied in the figures Id, 2d, 3d and 4d. In all 
the four figures it can be seen that the likelihood contours are parallel lines for all the three confidence levels. 
Unavailability of closed contours hinders our analysis considerably and due to this unbounded nature we are 
not able to put any finite bounds on both wo and cq from the present study and with the given choice of 
parameters. Looking at the figures, mathematically, one can obviously argue that the bound for the parameters 
are (— 00 , 00 ), but we do not think that it is physically acceptable. So, may be we have to devise some other 
methods to constrain the parameters for the logarithmic model, which we keep as an open issue for the time 
being. 

On a general note, it is observed that for Linear and CPL model, with the increase in the value of w 0 , the 
value of cq decreases. But in case of JBP model, an exactly opposite nature is witnessed. Increase in wo is 
accompanied by an increase in cq. It must also be noted that for both the Linear and CPL model w 0 lies in 
the negative range, whereas cq takes positive values. For JBP model, the scenario is a bit different. Both w 0 
and cq lies in the negative range. These results and the trends obtained for the involved parameters are quite 
different from the ones obtained by Biswas et al. in [68J using Loop quantum cosmology as background gravity. 
In [68], it was found that JBP model allows a wider range of values for the parameters compared to the linear 
and CPL models. Moreover with the addition of BAO peak analysis and CMB term the reduces the lower and 
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upper freedom of ui o and ui\. But in the present study no such trend is observed. 

In the fig. 5a, distance modulus /. i{z ) is plotted against redshift parameter z for all the four parametrization 
models and then compared with Union2 data Sample. The data sample is represented by a scatter diagram and 
the trajectories represent the n(z) vs z plots for various models. It is seen that barring the Logarithmic model, 
all the other models roughly fits the dataset. At lower redshifts (z < 0.3) it is seen that the fit is not quite as 
good as that at higher redshifts. The best fit is obtained around z = 0.4. Moreover for the Logarithmic model, 
there is almost no fit with the data. In fig. 5b, the trajectories are compared with SNe Type la 292 data. 
Here also it is seen that the logarithmic model does not fit the dataset whereas the other three models gives 
satisfactory fit with the data. In this case, the fit is considerably improved at lower redshifts compared to the 
union 2 data sample. At higher redshifts the curve fitting is proper. The best fit is obtained around 2 = 0.5. 
Although the best fit point has shifted to a higher redshift in case of SNe Type la 292 dataset compared to the 
Union2 dataset, yet the Sne data must be preferred compared to its counterpart because it gives a better fit 
with the parametrization models at lower redshifts with the given choice of parameters. Alberto et al. in |69] 
used Bayesian analysis to predict that JBP parametrization is preferred to CPL. But in this assignment, we did 
not find any result, which will either support or criticize the prediction. Finally, it must be stated that from the 
present study, we find that the predicted theoretical Logarithmic model does not fit either of the observational 
data. Moreover, for the logarithmic model, no finite bounds for the parameters could be obtained. 

In fig. 6, we have plotted the deceleration parameter q vs the redshift parameter z. Recently there have 
been quite a few works regarding cosmological deceleration acceleration transition mm In these works it can 
be seen that the transition from deceleration to acceleration occurs at z = 0.74 ± 0.05. We have obtained the 
expressions for the deceleration parameter for all the parametrization models using their best fit H ( z ) — z values 
of the parameters and then plotted them against redshift in fig. 6. From the figure it is seen that the transition 
from deceleration to acceleration regime occurs at z = 0.68, z = 0.74, 2 = 0.98 and z = 0.67 for the linear, 
CPL, JBP and logarithmic models respectively. So the result for the CPL model is in complete accordance with 
the works of Farooq et al. mum. Moreover the the values for the linear and logarithmic models almost lie in 
the redshift spectrum of transition obtained in Farooq et al. Only in case of the JBP model the transition is 
early which contradicts the result obtained by Farooq et al. Finally it must be pointed out that since we have 
fixed the free parameters with particular values, we have lost the overall generality of the study upto a certain 
extent. Correspondingly the EoS contours obtained could have been much broader, if these particular choices 
could have been avoided. So the present study might have lost some of its numerical accuracy. Nevertheless its 
cosmological significance is unquestionable. 
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